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Abstract 

A recent approach to the control of underactuated 
systems is to look for control laws which will induce 
some specified structure on the closed loop system. In 
this paper, we describe one matching condition and 
an approach for finding all control laws that fit the 
condition. After an analysis of the resulting control 
laws for linear systems, we present the results from 
an experiment on a nonlinear ball and beam system. 

1 Underactuated systems and 
the matching condition 

Over the past five years several researchers have pro- 
posed nonlinear control laws for which the closed 
loop system assumes some special form, see the con- 
trolled Lagrangian method of [g[ 3, |l0f the general- 
ized matching conditions of [l 2|, |13|, the inter- 
connection and damping assignment passivity based 
control of Q|, the A- method of || ||], and the ref- 
erences therein. In this paper we describe the im- 
plementation of the A-method of || on a ball and 
beam system. For the readers convenience we start 
with the statement of the main theorem on A-method 
matching control laws (Theorem 1). We also present 
an indicial derivation of the main equations. We 
then prove a new theorem showing that the family 
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of matching control laws of any linear time invariant 
system contains all linear state feedback control laws 
(Theorem 2). We next present the general solution 
of the matching equations for the Quanser ball and 
beam system. (Note, that this system is different 
from the system analyzed by Hamberg, fil|| .) As al- 
ways, the general solution contains several free func- 
tional parameters that may be used as tuning param- 
eters. We chose these arbitrary functions in order to 
have a fair comparison with the manifacturer's linear 
control law. Our laboratory tests confirm the pre- 
dicted stabilization. This was our first experimental 
test of the A-method. We later tested this method on 
an inverted pendulum cart, ||. 
Consider a system of the form 



grjX-i + [j k, r] iP x 



dV 

dx r 



(1) 



r = 1 , . . . , n, where gij denotes the mass- matrix, C r 
the dissipation, V the potential energy, [ij, k] the 
Christoffel symbol of the first kind, 



\jk,i] 



dgti _ dg jk 
dqi dq l 



(2) 



and u r is the applied actuation. To encode the fact 
that some degrees of freedom are unactuated, the ap- 
plied forces and/or torques are restricted to satisfy 
PJ-g^Uk = 0, where P? is a <?-orthogonal projection. 
The matching conditions come from this restriction 
together with the requirement that the closed loop 
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system takes the form Second, AJP^ and gij satisfy 

~ / gyl 

trjSP + [Jk^Wx k + C r + |K = o, r = 1, . . . , n, W + «] ^ - [rj, i]A* 



for some choice of g, C, and V. The matching condi- 

ti ° nSread , lp rflSU,. g(A^j , ~ g(A^ 

-V-n ^rr + 5 h ge m — j— — 



«\i \ 

-9ir-^ + [ij,r}K-[sjJ}4) =0, 



(6) 



y oq oq j Although the proof of this proposition may be found 

m |§j O' ana - 0' l° r convenience, we include an 
where r& is the Christoffel symbol of the second kind, indicial notation derivation of equations (^) and (7). 

Substitute equations (0) , (Q) for both and T^- into 
r*. = g ki [ij,£\. (4) tnc nrst °f equations (||) and multiply the result by 

the scalar 2 to obtain: 

If the matching conditions (0) hold, the control law prs. . fi-z„. Ftn.. 

wdl be given by ^k9 ~g q T ~ ^k9 ~g q l ~ ^9 -g^ 

u r = g rk {T k J -T k J )q l qi + [C r -C r ) - F k9 -g~r ~ F k 9 - P k 9 -r~ ■ 

(i-idV _ k {dv\ ^ Multiply by g r t and use that P is self-adjoint i.e., 

pr U d 9ij p r U <%hi pr x I d 9]t 

The motivation for this method is that H = t >' g q e t A r g^j r t A r g^i 

?9ij<l ll i J + V is a natural candidate for a Lyapunov g g g (8) 

, - _ . j-. _ pr <J 9v_ _ pr u 9ri _ p r u 9jr 

function because = —gijCq 1 . Following Q6|, in- ~ t <j r r t g^j r t g^ ■ 

troduce new variables = g^g^. We have 

_ _ USe ^* A r^T — 9^3 ff«t 9^3 ancl 

Theorem 1 The functions gij , V, and C satisfy (3) ^ 

m a neighborhood of xq if and only if X r gu = 9ri (9) 

, ^ v in (j§) to obtain (7). To derive (^), first, differentiate 

K [9^^ - g kl aj = 0, (g) with respect to qi to get 

pr( r M dV ~ki d v\_ n ^dg H _dg n dX e r 



Substitute equation ( |l0| ) into equation (||) and obtain 



and f/ie following conditions hold. First, there exists 
a hypersurface containing xo and transverse to each 

of the vectorfields X^P^dldx 1 on which q,-, is invert- v r ( ®\ , ^ ® A r . 
•L7 7 , . i j- ' -C F t + 9ij^r- + K 

ible and symmetric and satisfies \ oq J oq 1 



ible and symmetric and satisfies 

P, 



g ki Pl? = XlPjfgji. * yg q j + g q i g q j g q i + g q 



(11) 

dg r i . dg rj dg„ dg ]r dgij\ 
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Multiply by -P k Xl, use (g) and (0) to obtain 



(12) 



Finally, to obtain @, add to equation ( B ) an equa- 
tion obtained from ( |l2| ) by interchanging fe and t, r 
and s, £ and i. 

2 Matching and constant coef- 
ficient linear systems 

In this section, we prove that for linear time invariant 
systems any linear full state feedback control law is 
a solution to the matching equations. 

Theorem 2 When applied to linear, time- 
independent systems, the family of matching 
control laws contains all linear state feedback laws. 

Choose coordinates q % so that the desired equi- 
librium is at the origin, V = Vijq l q 3 + v k q k , and 
Ci = Crjq 3 , where g lj} Vij, v kl Qj, and P k are con- 
stant, and PJ1 has rank n u . Clearly, V — Vijq l q 3 , 
and Ci = Cijq 3 is a solution to the matching equa- 
tions when gij, Vij, and Cij are constant provided 
g-ij and Vij are symmetric, P£ (^g kt Vij — g kl Vij^j = 0, 

and PI (g ki Cij - g kl C l ^j = 0. Let u k = v k + a kz q l + 

b k iq l be an arbitrary linear control law, satisfying 
PJlg kl ui — 0. Comparison with equation (3) gives 



and 



Thus, 



and 



9r k (V'i;, ;f\,. 



9rk ( K 9 kl C'ij - g hl Cij^ 



Vij = gi P g pr {V r] - a rj ) 



C/j — gi P g pr (C r 



It remains to check that we can find a symmetric, 
nondegenerate matrix ~g k% so that the resulting Vij is 



also symmetric. The symmetry of Vij will follow if 
we have 

ge P g pr (Vr] - a rj ) - dj P g pr (V r i - a r e) = 0, 

and, therefore, we need to find a symmetric, nonde- 
generate matrix gi p satisfying this equation. The ex- 
istence of such matrix is guaranteed by the following 
simple observation. 

Lemma 1 Given any real n x n matrix R, there is 
a nondegenerate symmetric matrix X so that 

RX - X T R T = 0. 

Indeed, setting X — QYQ T , results in the following 
equation for Y: 

Q~ l RQY - Y T {Q- 1 RQ) T = 0. 

Hence, without loss of generality we may assume that 
Q~ 1 RQ is a real Jordan block (see pi). 

f X 1 . 
A 1 . 



or 



(a 
b 
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-b 1 
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In each case Y 



/0 




1 \ 







solves the 



yi o ... o J 

equation. 

Note that the result of Lemma 1 is true for matrices 
with coefficients in any field. This is proved in Jig] , 



3 Example: 
Beam 



The Ball and 



In order to demonstrate the approach described 
above, we have implemented one of the control laws 
from the family of control laws described in the first 
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Table 1: The physical parameters of the system. 



4 = 


length of the beam 


= 0.43m 


Ib 


= |msrg = ball inertia 


= 4.25 x 10~ 6 Kg m 2 


h = 


length of the link 


= 0.11m 


h 


= inertia of the beam 


= .OOlKg m 2 


r 9 = 


radius of the gear 


= 0.03m 


Is 


= effective servo inertia 


= 0.002Kg m 2 


TB = 


= radius of the ball 


= 0.01m 


9 = 


= gravitational acceleration 


= 9.8m/s 2 


TUB 


= mass of the ball 


= 0.07Kg 


so 


= desired equilibrium position 


= 0.22m 


m>b = 


= mass of the beam 


= 0.15Kg 


CO 


= inherent servo dissipation 


= 9.33 x 10- 10 Kg m 2 /s 


mi = 


= mass of the link 


= QMKg 










Figure 1: The ball and beam system 

section on a ball and beam system, Figure 1 (this 
system is commercially available from Quanser Con- 
sulting, Ontario, Canada). 

The s— coordinate is unactuated, the 9— coordinate 
is actuated by the servo, and the objective is to bring 
the ball to the center of the beam. The physical pa- 
rameters of the system are given in Table 1. 
One can express a as a concrete function of 9 from 
the kinematic relation 



(4(l-cos(a))-r s (l-cos(0))r 
+ (£ b sin(a)+£i-r g sm(0)f 



The kinetic energy of the system is 



1 2 2 1 . 
T = -m b s a + -Ib{oi 

The potential energy is 

V = l; m i 9 r g sin (#) H 
+ tub g s sin(a) , 



1 



-ha 1 



-(m b +mi) g 4 sin(a) 



and the dissipation is C\ = 0, C2 = cq9 . After rescal- 
ing, we get 

(l-cos(a)-a 2 (l-cos((9))) 2 

+ (sin(a) + ai — a 2 sin(0)) 2 = a 2 , 

T= ^s 2 + ^(a 4 +(a 3 +5/2 S 2 )(a , (e)) 2 )e 2 +a'(9) s9 ,C, =0, 

V = a 5 sin(6») + (s + a 6 ) sin(a(6»)), 
C\ = 0, and Ci = 0,76, where the at are the dimen- 



sionless parameters, 



ai 



Li 

4' 

a 5 



a-2 



4' 

mir g 
2m B r B ' 



a 7 



a-3 



«(> 



(h + Ib) 
Ib 

_ £b(mb + mi) 
2m B r B 

Co 
m B 
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II 

Ib' 



The notation ' is used to denote a derivative of a 
function of one variable. For general underactuated 
systems, the use of the powerful A-mcthod to solve 
the matching equations is discussed in || ||]. For 
systems with two degrees of freedom, the A-mcthod 
produces the general solution to the matching equa- 
tions in an explicit form, When applied to the 
ball and beam system, the explicit family of control 
laws is given by equation (|J) with the following ex- 
pressions for 5, V , and C, where 



g n =^ 2 (a) %(*,*)) + 10 



512 



AWWW) 



-{911 - <rgu), 922 = -{912 - 0-512), 



C-2 
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Jo Viwmv) Jo 

C 1 = (g li 9 ll r 1 (c 1 -g 1 ^ 2 d 2 ), 

where 

MM) -2&8&, 

5s5i2 

<t(s,0) = /tii (a) - — Ai'iCa), 

OS 

y = ^(a)s - s + / ip(r) dr, ) dn). 

Jo LhW 

Here h(y), w(y), Hy(a) are arbitrary functions of one 
variable, and C 2 is an arbitrary function which is odd 
in velocities. 



4 Experimental Results 

Our experiments were conducted on the Quanser ball 
and beam system. The control signal is a voltage sup- 
plied to the servo and the sensed output of the system 
is s and sampled at 300 Hz. A Quanser MULTIQ® 
data acquisition card is used for the analog signal 
input and output. The velocities are computed via 
numerical differentiation using the forward difference 
algorithm. The control law produces a voltage signal 
and is supplied through the D/A converter to the DC 
servomotor via an amplifier. The relation between 
the control voltage, Vi n , and the torque, u (= u 2 in 
equation (|j)), is K^N^O, where R m = armature re- 
sistance = 2.6 17, N g — gear ratio = 70.5, K m 
motor torque constant = 0.00767 Volt-sec. 

Any stabilizing linear control law for this system 
is specified by four constants. The nonlinear control 
laws in our family are specified by the four arbitrary 
functions: fiy(a), h(y), w(y), and C2{s,9, s,0). We 



chose 

m(a) = 1.0849 exp(4.7845sin(a)) 
h(y) = 1.1031, w(y) = 0.0023y 2 , 
62(8,0,8,6) = -g u ■ (1 + s 2 + 1O0 2 )(-/xs + <t0). 

These functions produce the control law, u, in 
rescaled units. The values of the constants ay through 
a7 are as follows 



ay = 


0.2547 


a 5 = 


0.1889 


a.2 = 


0.0588 


a 6 = 


42 


03 = 


236.294 


a 7 = 


5 x 10 


04 = 


471.126 







The final control signal is obtained by converting 
back into MKS units and using the formula in the 
preceding paragraph to get the input voltage. These 
choices were made from the following considerations. 
The form of the function \xy was chosen to simplify 
the integrals in the expressions for y, ib, and gyy. The 
form of C2 was chosen to ensure that Cys+C20 would 
be positive (for H to be a Lyapunov function). Fi- 
nally, the coefficients in these functions were chosen 
so that the linearization of the nonlinear control law 
would agree with the linear control law provided by 
the manufacturer. 

Extensive numerical simulations done using 
Matlab® confirm that the nonlinear control law sta- 
bilizes the system. The linear control law appears to 
stabilize the system for a wider range of initial condi- 
tions than the nonlinear control law. This is an em- 
pirical observation, not a mathematical fact. Find- 
ing an adequate mathematical framework to compare 
different control laws is a very interesting unresolved 
problem, see |]. Usually, given two locally stabilizing 
control laws, there exist initial conditions stabilized 
by one but not by the other. For example, one set of 
physically unrealistic initial conditions with a large 
angular velocity 9 = 3.6 (or 158 rad/sec in physical 
units) is stabilized by our nonlinear control law but 
not by the linear one. 

We have implemented the nonlinear control law in 
the laboratory. The laboratory tests confirm the pre- 
dicted behavior of the nonlinear controller. Figures |^ 
and H show a comparison of the time histories of the 
ball position (s) and angular displacement (0) for the 
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-- Nonlinear control law 
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Figure 2: Ball position response 

linear and nonlinear control laws. In both cases the 
control signal reached the saturation limit for a short 
duration during the initial rise of the response. The 
difference in the steady-state values of the responses 
is attributed to a lack of sensitivity of the resistive 
strip used to measure the ball position. 

5 Conclusions 

The A-method produces explicit infinite-dimensional 
families of control laws and simultaneously provides 
a natural candidate for a Lyapunov function. When 
this method is applied to linear time-invariant sys- 
tems, the resulting family contains all linear state 
feedback control laws (Proposition 2). In this paper 
we also present the results of the first implementation 
of a A-method matching control law on a concrete 
physical device, the ball and beam system. The ex- 
perimental results agree with theoretical predictions 
and numerical simulations. In our experiments we 
observe that the linear control law performs better 
than our nonlinear control law for the ball and beam 
system. However, in a later experiment with an in- 
verted pendulum cart, S, we found that a properly 
tuned A-method matching control law performed bet- 
ter than the corresponding linear one. At the mo- 
ment, it is not known for which systems matching 
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Figure 3: Angular displacement response 

control laws will perform better. This is an impor- 
tant problem that must be resolved. 
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